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Abstract 

In this paper we present a solution to the /i problem in an SO(IO) super- 
Oh. symmetric grand unified model with gauge mediated and D-term supersym- 

r-| ! metry breaking. A Peccei-Quinn symmetry is broken at the messenger scale 

M ~ 10^^ GeV and enables the generation of the /x term. The boundary 



X 



conditions defined at M lead to a phenomenologically acceptable version of 
the minimal supersymmetric standard model with novel particle phenomenol- 
ogy. Either the gluino or the gravitino is the lightest supersymmetric particle 
(LSP). If the gravitino is the LSP, then the gluino is the next-to-LSP (NLSP) 
with a lifetime on the order of one month or longer. In either case this heavy 
gluino, with mass in the range 25 - 35 GeV, can be treated as a stable particle 
with respect to experiments at high energy accelerators. Given the extensive 
phenomenological constraints we show that the model can only survive in 
a narrow region of parameter space resulting in a light neutral Higgs with 
mass ~ 86 — 91 GeV and tan/3 ~ 9 — 14. In addition the lightest stop and 
neutralino have mass ~ 100 — 122 GeV and ~ 50 — 72 GeV, respectively. 
Thus the model will soon be tested. Finally, the invisible axion resulting 
from PQ symmetry breaking is a cold dark matter candidate. 



1 Introduction 

Supersymmetry (SUSY) is a strongly motivated candidate for new physics 
beyond the Standard Model (SM). It provides a natural framework for re- 
solving the hierarchy problem. The minimal supersymmetric standard model 
(MSSM), with a conserved R-parity, has an economical particle content with 
well defined interactions most of which are already constrained by experi- 
ment. It has two Higgs doublets {Hu and i/^); necessary for giving mass to 
both up and down quarks, respectively. In the MSSM, electroweak symmetry 
breaking (EWSB) occurs naturally since rn?H^, the soft SUSY breaking (SSB) 
mass of Hui is automatically driven negative as a result of a large top quark 
Yukawa coupling. 

The MSSM solves the hierarchy problem by allowing for dimensionful 
SSB parameters of order the electroweak scale and protecting scalar masses 
from large radiative corrections above the SSB scale. However, this in itself 
is not sufficient to solve the hierarchy problem. In addition, it is necessary 
to demand that the /x parameter (where /i is the bilinear Higgs coupling in 
the superpotential of the form ixHy^Hd) is also of order the electroweak scale. 
Consider the vacuum conditions obtained by minimizing the tree level Higgs 
potential, we have 

^2_ M| , ml^- ml Jan? (5 

where tan/3 = {Hy) / {Ha) , sin2(3 = 2B/m\, niA is the mass of the CP odd 
Higgs and B is the SSB Higgs bilinear coupling. On the left hand side of Eq. 
m, the /i parameter, multiplying a supersymmetric fi term in the Lagrangian, 
breaks no SM symmetries; it could in principle be as large as the Planck or 
GUT scales. On the right side, the Z boson mass and the SSB Higgs masses 
are of order the electroweak scale. Clearly, all three scales, i.e. /i, Mz and 
the SSB scales must be of the same order. This is the n problem ^. 

In order to avoid large values for /i, a symmetry is needed which prevents 
the n term at tree level, but allows such a term once this symmetry is broken. 
Moreover, since in the MSSM the fi term is contained in the superspace 
potential, there are two possibilities - 1) it can be generated via a term in the 
Kahler potential (at tree level or radiatively) once supersymmetry is broken 
or 2) no supersymmetry breaking is required if it is generated through higher 



dimension operators in the superspace potential. Several simple mechanisms 
for generating a /i term have been suggested 0, Q . 

In the context of gauge mediated supersymmetry breaking (GMSB) mod- 
els [^, there is an additional problem, fi can be generated at one loop order 
once supersymmetry is broken. However, the B parameter (the SSB scalar 
Higgs bilinear coupling) is usually generated at the same loop order and is 
too large. A solution generating B at higher loop order than /x was proposed 
ini. 

In this paper we use an extension of the GMSB model discussed in Ref. p 
to solve the /i problem. This model has GMSB with Higgs-messenger mixing 
in an SO(IO) theory and naturally leads to a gluino LSP. The gluino LSP is 
stable due to R-parity conservation. The specific signature of a gluino LSP 
i.e. missing momentum has been analysed in Ref. 0] for LEP and CDF data 
and in Ref. |^ for CDF data. Ref. |^ concludes that a stable gluino with 
mass in the range 25 — 35 GeV is still allowed by both the LEP and CDF 
data. Our model, with a modest adjustment of parameters, gives a gluino 
with mass in this range. 

The /i term is absent in this model (at tree level) due to a U(l) Peccei- 
Quinn (PQ) symmetry. Both SUSY and PQ symmetry are broken when the 
chiral superfield X develops a vacuum expectation value (vev) 

{X) = M + e^Fx. (2) 

Note, the following dimension five operator in the Kahler potential 

K D -^XhOJ, + h.c. (3) 

is however allowed by the symmetries. Thus we find /i ~ Fx/Mp. B, on the 
other hand, is generated radiatively via renormalization group (RG) running 
below the messenger scale M.Q 

In addition to solving the fi problem, the PQ symmetry provides a natural 
solution to the strong CP problem 0. The strong CP violating 6 term 
dynamically tracks to zero. Moreover, as a bonus, the axion is a candidate 
for cold dark matter. 



^Note, the niechaiiisin used here to generate the /z term is a combination of the ideas 
discussed in 0, 01 • 



In the next section we discuss the model, saving some of the details for the 
appendices. We derive the low energy spectrum consistent with electroweak 
symmetry breaking, gauge coupling unification and third generation quark 
and lepton masses. We then consider experimental constraints which con- 
strain the available parameter space to a very narrow region. In this region 
we find a light neutral Higgs with mass ~ 86 — 91 GeV and tan/3 ~ 9 — 14. 
In addition the lightest stop and neutralino have mass ~ 100 — 122 GeV 
and ~ 50 — 72 GeV, respectively. Finally, in an appendix we investigate the 
possibility of obtaining a reasonable mass for the tau neutrino in the model. 
Clearly this model is preeminently testable. 

2 The model 

The theory at the GUT scale is defined by the SO(IO) invariant superpo- 
tential W D Wi + W2 + W^ and a non-renormalizable term in the Kahler 
potential K where 

Wi = 16310^163, 

W2 = XalOnAlOA + AxXlOi, (4) 

W3 = Xif]iAr]i + \2f}2Arj2 + XXf]ir]2. 

K DA;,-— 10^ + /i.e. (5) 
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(I63, rji, r]2) are 16's, (r^i, 1)2) are 16's, {'^Oh, 10^) are lO's, (X) is a singlet 
and {A) is an adjoint under SO(IO). 

At the GUT scale, the theory is invariant under a U(l) PQ and an R 
symmetry. The R symmetry is broken spontaneously at the GUT scale. The 
PQ symmetry, however, is not broken at the GUT scale and prevents a ix 
term in the superpotential. The PQ and R charges of the fields are defined 
in Appendix ^. 

Wi contains the coupling of the third family matter multiplet (I63) to 
the Higgs field (10//) which includes both the weak doublet and color triplet 
Higgs fields. 



W2 serves two purposes. In the first case, it provides doublet-triplet 
splitting using the Dimopoulos-Wilczek mechanism [T^. The adjoint field A 
gets a vev 

{A) = {B- L)Mg, (6) 

where B — L (baryon number minus lepton number) is non- vanishing on color 
triplets and zero on weak doublets and the singlet X gets a vev 

{X) = M + d''Fx. (7) 

This gives mass of order Mq to the color triplet Higgs states and of order M 
to the weak doublets in 10a- The Higgs doublets in 10^ remain massless. 
The SUSY breaking vev, on the other hand, exhibits the second purpose for 
W2. 

In the second case, W2 and PV3 also provide the messengers for SUSY 
breaking.0 The auxiliary field 10^ and the fields f]i.,r]2 feel SUSY breaking 
at tree level due to the vev Fx- They are thus the messengers for GMSB 
[|, [ll|. We take the messenger scale M ~ 10^^ GeV with the effective SUSY 



breaking scale in the observable sector given by 

K = Fx/M ^ 10^ GeV. (8) 

Thus the Higgs field in this model plays a central role with regards to 
supersymmetry breaking. It is this central role which also provides a natural 
framework for solving the /i problem using the PQ symmetry. When X 
gets a vev, both SUSY and the PQ symmetry are broken. The /x term is 
generated at the scale M 

. = A.^. (9) 

while B remains zero at tree level. 

The PQ symmetry solves the strong CP problem and produces an axion; 
the Goldstone boson of the broken PQ symmetry. The axion gets mass 
due to the QCD chiral anomaly of order m\ = ^mlN'^ (i+z)^ EH where 



Z = mu/md ~ 0.56, fa = M = 10 GeV is the PQ symmetry breaking 

^Due to an accidental cancellation, gluinos receive no mass at one loop from W2 ■ Thus 
W3 is introduced with additional messenger fields (??i, ?yi, ?72, ?y2) contributing to the masses 
of gauginos and scalars at the scale Mq- 



scale and A^ = 3 is the number of families. Putting in the numbers we find 
ma ~ 2 X 10-^ eV.0 

We refer the interested reader to Appendix ^ for the complete model 
defined at the GUT scale. Note, in order to obtain realistic t, b and r 
masses, we find it necessary to abandon Yukawa unification at Mq- How 
this is obtained in the complete model is discussed in Appendix |B[ Finally, 
a simple extension of the model to include a r neutrino mass is presented in 
Appendix 0. 

3 Boundary conditions at the messenger scale 

The boundary conditions at the messenger scale are determined by two 
sources of SUSY breaking, gauge mediation and D-term 0. The messen- 
gers give mass to the gauginos and Higgs at one loop and to squarks and 
sleptons at two loops. Since the color triplet messengers have mass of or- 
der the GUT scale, the gluino mass is suppressed compared to the other 
gauginos. The gaugino masses (at M) are given by 



mg = f A62, 

M2 = fA(l + f62), (10) 

Mi = |^A(l + 462). 
where 

The two loop GMSB contribution to the scalar masses is given by 

(12) 
where a = ^^^ , C3 = 4/3 for color triplets and zero otherwise, C2 = 3/4 
for weak doublets and zero otherwise and Ci = (|)(F/2)^. a and b are free 
independent parameters which we use to fit the data at the EWSB scale. 



■^We take M = 10^^ GeV, so that the energy density in the axion field does not over 
close the universe. 



The gravitino mass is given by 

where A^usy is the scale of SUSY breaking and Mp = (SttGw)"^/^ = 2.4 x 10^^ 
GeV is the reduced Planck mass. For A'j.^^y = Fx — 10^^ GeV^ the gravitino 
mass is 

ma = ^ A - 0.024 GeV, (14) 

making the gravitino the LSP. However, this conclusion is model dependent. 
For example, we show that in the particular model of SUSY breaking 
discussed in Ref. [|^, the gravitino mass is significantly larger. In this 
model the field which gets both a scalar and F-component vev is the third 
component of an SU{2)p vector field 5*3. In this theory X is a composite 
field with (X) = M = Sl/Mst and Fx = 2S^FsjMst. The gravitino mass is 
therefore given by |[Tl[ 



^G = -JkVT- = TTT^ — T7 ^' (15) 



Fs, ^ 1 VMM Tt 
V^Mp 2VS Mp 

with A still given by A = Fx/M. The gravitino mass is thus enhanced by the 
factor \J ^j^- For example, letting the string scale Mgt = Mp and requiring 
the scale of PQ symmetry breaking M = 10^^ GeV, we find rriQ = 18.6 
GeV.0 To conclude, in this model, with A = 10^ GeV, the gravitino, gluino 
and wino masses (at M) are given by 

m^ = 18.6 GeV 
mg = (A)2 X 14 GeV (17) 

Mo = 340 GeV. 



^It is necessary to check that the supergravity contribution to squark and slepton 
masses is small compared to the GMSB contribution. This ratio scales as 



M2 V3a2 Mp ■ 



(16) 



Thus S3 cannot be much larger than 10^^ giving ^^ ~ 0.04. Taking the string scale 



Ah 
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Alst — Mp and requiring the scale of PQ symmetry breaking M = 10 GeV, one gets 
53 ~ 1.5 X 10^^ and ^^ ~ 0.056 which is reasonable. 



Hence, either the gluino or the gravitino is the LSP depending on the par- 
ticular SUSY breaking model and the value of the parameter 6.0 

For phenomenological reasons we assume that SUSY is also broken by 
the D-term of an anomalous U{l)x gauge symmetry as already discussed 
in Ref. 0. Moreover, the GMSB and D-term contributions are necessarily 
comparable.0 The D-term contribution to scalar masses is given by 



6DK = dQ^Mi, (18) 

where Q^ is the U{l)x charge of the field a and d is an arbitrary parameter of 
order 1 which measures the strength of D-term versus gauge-mediated SUSY 
breaking. The value of Q^ for a = 16, 10, 1 of SO(IO) is given by 1, —2, 4 

i- 

We now summarize the messenger scale boundary conditions. The gaug- 
ino masses are given by Eq. [1^ and scalar mass by Eqs. |12| and |18|.[] D-term 



SUSY breaking only contributes to the scalar masses. The SSB trilinear 
scalar coupling A and scalar Higgs mass squared B vanish at tree level but 
are generated via RGE running below M. We have chosen a = 10~^. We 
then determine the free parameters A, b, d and /i (at M) and acuT, ^3, 
the top (At), bottom (A^) and r (A,-) Yukawa couplings (at Mq) by fitting 
the low energy data which we take to include mj, m^, m^, aem, «« and 
sin^ 6*1^-0 Imposing gauge coupling unification at the GUT scale, we renor- 
malize the effective Lagrangian parameters to the EWSB scale using one 
(two) loop equations for dimensionful (dimensionless) parameters. We have 
also included the effect of a light gluino in the running of rrib and as below 



the EWSB scale [O]. Finally the one loop SUSY threshold corrections for 



the b and r masses at the EWSB scale have also been included. 



^In order to have b = 0.1 with M/Mq -- 10"'' we need to take A2/A1A2 ~ 10^. 

^ Ref. |1^ presents a model which dynamically breaks SUSY and leads to comparable 
SUSY breaking effects from gauge-mediated and D-term sources. 

^Note, the Higgs states also get a SSB mass correction at one loop due to their direct 
interaction with the messengers. This correction is negligible however because of the small 
values we have chosen for the free parameters a, b. 

^ Note, we allow for a small one loop threshold correction to gauge coupling unification 
at Mg which we have parametrized by the free parameter 63 = (a^ — aG)lotQ evaluated 
at Mg- 



4 Phenomenology at the EWSB scale 

The parameter b is varied to keep the gluino pole mass at 30 GeV. | d and A 
set the scale for squark, slepton and gaugino masses. We have examined cases 
of fixed A, varying d and vice versa in order to study the effects of these two 
SUSY breaking mechanisms on the low energy phenomenology separately. 
tan/? is solved from Eq. ^and helps determine the quark and lepton masses. 
We have allowed for values of (63 < 4%). 

In the first part of our analysis, A is fixed to the value A = 10^ GeV while 
d is varied. Fig. |l| shows the values of |/i| at the messenger scale giving the 
best low energy fit. Note that fi < and |yu| increases with d. The reason 
for negative n can best be seen by considering the following equation 

At small tan/3, the second term on the right hand side of the equation is 
negligible. Although rn^j^ and m\j have a common value at the messenger 
scale, Tn\j^ is always driven to smaller values than m\j by the RGE because 
of the larger top Yukawa coupling. We also know that tan 2/3 < as long 
as tan/3 > 1. We therefore need i? > at the EWSB scale. The one loop 
/3-function for B is given by 

/3i^) = B{m? + 3|nP + \Yr? - 3^i - lal) 

(20) 
+/i{6Ar/ + QA,Y^ + 2AX + QglM^ + |(7?Mi}. 

Since B and the trilinear couplings At^b.r vanish at the messenger scale, we 
must choose yU < in order to get a positive B at the EWSB scale. 

The values of |/i|,A/B and |-B//i| at the Z scale are also shown in the 
plot. It is notable that RGE running gives small values of 5, eliminating 
fine tuning in the Higgs potential, and thus giving a good solution to the /x 
problem. Note, |/i|, vB are both increasing functions of d] with a sharper 
rise at small values of d. This can be understood using Eqn. |l[ For moderate 

^Using the RGEs we evaluate m^^^ at rrig and then calculate the one loop corrected 
gluino pole mass. It is the pole mass which is constrained to lie between 25 and 35 GeV. 
Note, that an 18 GeV MS running mass defined at Mz is roughly equivalent to a 30 GeV 
pole mass. 



values of d, the M| and rn^j terms are negligible (the latter due to the fact 
1 771^ I < |m|f I and the factor of tan/5^ — 1 in the denominator). Hence the 
approximate relation /i^ = —m?^ > holds. Increasing d leads to a linear 
increase in \m'jj \ (the U{l)x charge of H^ is -2) and therefore in |/i|. For 
very small values of d, \rn?fj \ becomes comparable to M|/2 and the significant 
cancellation in the relation /x^ k. — M|/2 — m?H^ leads to a sharp decrease 
in /i. Since B is generated from /x via RGEs, the dependence of a/B on d 
follows that of /i. 

We mentioned that fn^j^ — mfj is zero at the messenger scale but is 
negative at the Z scale. Decreasing d has a small effect on the value of 



rriH^ — iTijj at the Z scale. However from Eq. |T9| we see that a sharp 
decrease in B at small values of d has to be compensated by small values of 
I tan2/5|. This is the reason for the sharp increase in tan/?, plotted in Fig. ^ 
at low values of d^ 

In Fig. ^ we plot the masses of the Higgs states. The masses of A, H^ and 
if"*" are determined at tree level while we have included the one loop SUSY 
threshold corrections to the mass of h. Since B increases with d, the masses 
of A, H^ and H^ also increase while A stays the lightest. One might think 
that since the mass of the lightest Higgs state h is an increasing function 
of tan/?, it should decrease when increasing d. However since our model is 
constrained and B increases with d, the effect of B dominates over that of 
tan/? and the mass of h slowly increases with d. The reason for the slow 
increase is that the mass of h is determined by EWSB and hence it cannot 
be a strong function of d. 

In Fig. ^, we have magnified the small d and A regions in order to show 
the rapid variation in the value of nth here. The solid line shows the change 
in the mass of h versus d for a fixed A = 10^ GeV while the dashed line shows 
the variation versus A with a fixed d. As discussed above, B decreases with 
decreasing d. This decrease in B is directly reflected in the decrease in the 
tree level value of nih- Similarly the value of B decreases when A decreases 
and thus so does m/i.|^ 

^°Note, a small | tan2/3| gives a large tan/? assuming 7r/4 < f3 < 7r/2. 

^^The Higgs mass apparently increases without bound as A increases. We believe this 
result is due to the fact that we only use the tree level Higgs potential for EWSB. Note, 
we have checked that we cannot obtain reasonable fits to the data for A > 1.1 x 10^ GeV; 
hence we believe that our range for the allowed Higgs mass will not be significantly affected 
when including higher order corrections to the Higgs potential. 



For completeness, we show the behavior of tan/? with changing A and 
fixed d = 1 in Fig. ^; tan/3 decreases rapidly for A < 10^ GeV. This is because 
I tan2/5| must increase to compensate for |m^ — m'jj \ which is decreasing 
proportional to A^ (see Eq. |19]). 

Fig. ^ contains plots of the masses of charginos and neutralinos versus 
d. The mass of the heavy chargino [xt) ^^^ the two heaviest neutralinos 
(Xg, X4) increase, as |//| increases with d. However, for most values of d, the 
mass of Xi scales as Mi while the masses of X2 ^^^ xt scale as M2 and do 
not run with d. At very small d, tan/3 gets very large and the off diagonal 
elements, proportional to sin (3 in the chargino and neutralino mass matrices, 
become larger than the diagonal elements including fi resulting in a sharp 
drop in the masses. 

Squark and Slepton masses are plotted in Fig. |^. The D-term contributes 
positively to the masses of squarks and sleptons since the U{l)x charges of 
these fields are +1. Consequently, we see an increase in their masses with 
increasing d. The mixing due to A and fi terms for the third generation 
squarks and sleptons is very small. Nevertheless, due to the different bound- 
ary conditions at M, the right handed squarks and sleptons are lighter than 
the left handed ones. The right handed stop is always the lightest and below 
d ~ 0.7 it becomes lighter than the top. This is because the right handed 
stop mass squared is driven negative by RGE as a consequence of the large 
top Yukawa coupling. 

Finally in Fig. |^ we plot the masses of squarks, sleptons, charginos, neu- 
tralinos and Higgs fields versus A. As one expects, all these masses increase 
with A except the lightest Higgs whose mass is determined by EWSB. Re- 
call that all scalar and gaugino masses at the messenger scale are directly 
proportional to A. 



5 Laboratory and cosmological constraints 

Let us first consider the heavy gluino. If it is the LSP it was shown that it 
can survive in a narrow window with mass between 25 and 35 GeV |^, |^ . We 
note that to get this limit, Refs. 0, H assume very large squark masses. Ref. 
[§] however also argues that lowering the squark masses increases the allowed 
range for the gluino mass. Now consider the possibility that the gravitino is 
the LSP and the gluino is the NLSP. In this case we must check whether the 
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analysis of Refs. [0, |^ still applies, i.e. whether the gluino lifetime is greater 
than ~ 10~^ s. The decay rate of the gluino to a gluon and a gravitino is 
given by 

r "' ^i (1 <f (21) 

Even in the case that A^^^^ = Fx — 10^^ GeV^, which gives a very light 
gravitino as in Eq. |T^, a gluino mass of 30 GeV has a lifetime of Tg ~ 
2 X lO^s ~ 1 month. We therefore conclude that the gluino is in all cases 
a stable particle with regards to detector experiments. Hence the missing 
momentum analysis of Refs. [^ |^ is relevant and a heavy gluino NLSP with 
mass in the range 25 — 35 GeV is still viable. 

Note there are several significant advantages for having a heavy gluino 
LSP or NLSR 

• It reduces the fine tuning necessary for EWSB, since the dominant 
contribution to scalar masses due to RG running from Mq to Mz comes 
from color corrections proportional to the gluino mass squared p!5[1. 



Even if the gluino is the NLSP, its lifetime is long enough for it to be 
a candidate for the UHECRon, i.e. the source of the ultra high energy 
cosmic rays p6|] . 

The model with a Higgs mass of order 90 GeV and a stop mass less 
than the top satisfies some of the dynamical constraints necessary for 
electroweak baryogenesis in supersymmetric theories |T^ . 



We now consider the LEP constraints on other SUSY parameters in our 
model. The most important constraints come from the latest Higgs search 
results at LEP [|l^]. The light neutral Higgs h and the CP odd Higgs A in 
the MSSM are produced at LEP via the Higgs-strahlung process e~^e~ —>■ hZ 
or the pair production process e~^e~ — *■ hA. h and A decay predominantly 
into bb and t^t~ . Thus LEP experiments search for either a hb or t^t~ 
plus the decay products of the Z; or for bbbb and r^r^bb. In our model the 
off-diagonal elements in the stop mass-squared matrix are very small, thus 
the LEP limits for the neutral Higgs in the no stop-quark mixing scenario 



are most applicable |T^, |T9[. These limits are very severe. Looking at our 



data in Figs. and ^, it can be seen that only d ~ 0.40 — 0.45 survives the 
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LEP constraint for A = 10^ GeV. With these values of the parameters, the 
mass of the hghtest neutral Higgs resides in the narrow range ~ 86 — 91 GeV 
with tan/3 ~ 9 — 14. At this point, our model also survives the limit on the 
mass of A as indicated in Ref . [|1^ . H~^ and H^ are also too massive to be 
constrained. For d ~ 0.40 — 0.45 and A = 10^ GeV we find the lightest 
stop and neutralino with mass in the range 100 - 122 GeV and 50 - 72 GeV, 
respectively. Note, if A increases, then the Higgs mass and tan/? remain 
unchanged while all other masses increase .Q 

This narrow region of parameter space is obtained with the fit values of 
the parameters aa = 5.28 - 5.49 x 10^^ GeV, eg = 2.44 - 2.42%, A^,*,^ = 
(0.065, 0.42, 0.096) - (0.042, 0.42, 0.061) at Mg and /i = (-85.4) - (-139) 
GeV and b = 0.06 at M. Note, we are not able to get a good fit to the data 
assuming Yukawa coupling unification. This was not the case in Ref. and 
is due to the fact that, unlike the model presented here, the parameter B 
was a free parameter. In the appendix, we show how to complete the model 
described earlier in order to accomodate non-unification of Yukawa couplings 
in this SO(IO) theory. Clearly this feature of the model is not very satisfying. 

The process e+ e~ —>■ hadrons can constrain the chargino mass based on 



the OPAL 2 a bound on new physics at ^/s=172 GeV [21 . We analysed the 
contributions to e"^ e~ -^ hadrons coming from chargino and neutralino pair 
production using {SPYTHIA, A Supersymmetric Extension of PYTHIA 5.7 
[^), which has been modified by S. Mrenna and K. Tobe to accommodate 
a gluino LSP. The result is that except for very small values of d or A which 
are strongly ruled out by the Higgs constraint, the cross section is below the 
limit of 3 pb reported in Ref. |^ .0 

We have evaluated the rate for b -^ sj.f^ The ratio of the SUSY am- 
plitude to the SM amplitude for (varying rf-fixed A) and (varying A-fixed d) 
are given in Fig. 0. The main SUSY contribution in our model comes from 



^^Note, a recent LEP bound on a heavy gluino LSP using stop production and decay 
p(| does not constrain the model since the stop mass in our case is larger than the values 
probed in this search. 

^•^For example, when A = 10^ GeV, the cross section of the process e"*" e~ -^ hadrons 
through charginos and neutralinos exceeds the 3 pb limit only for d < 0.37. 

'^^In order to do this calculation we need to assume some values for flavor mixing ele- 
ments. As a rough estimate we use the observed CKM matrix elements in the appropriate 
places. This is clearly just a rough estimate which gives us at best an order of magnitude 
approximation. 
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the charged Higgs-top loop. Fig. ^ shows that we always get an amplitude 
larger than the SM amplitude approaching it at large A since /i < 0. Note, 
our result is at the limiting edge of values allowed by CLEO data P3 . Our 



result is however incomplete. We approximate flavor mixing using the known 
CKM elements as a crude approximation to squark - quark flavor mixing and 
we only use a one-loop analysis. With regards to the latter, it is clear that 
a one loop analysis is not sufficient since the higher order corrections have 
recently been shown to be very important (see for example Ref. p^)- Hence 



within the approximations considered here, the model is roughly consistent 
with the observed rate for 6 — > S7. 

Now consider the low energy consequences of the PQ symmetry which 
is spontaneously broken at the scale /a = M = 10^^ GeV. This generates 
an invisible axion which is a significant component of the energy density of 
the universe and thus a candidate for cold dark matter [^. The axion is 
predominantly the angular part of the complex scalar field component of 
the chiral superfield X. The radial part of this scalar field and the fermion 
component of the supermultiplet are named saxino {S) and axino {A) in 
the hterature, respectively. In our model, these fields obtain large radiative 
masses via a messenger lO^i loop. 

•^'A = ^A. Ml = ^A^ (22) 

The axino decays mainly to a gluino and a gluon. The decay rate is 
calculated to be 1261, 



r,^,, = ^A.^)^. /(x) = -^ + ^, (23) 

and we find IQ-'^'^GeV < V .- < IQ-^^GeV for 0.2 < ^^^ < 1 and 0.5 < 

Xx < 3. This translates into an axino lifetime of 10~^s < Ta^~ < 10~^s. 
The saxino decays either to 2 gluons or 2 gluinos. The decay rate of the 



saxino to 2 gluons is given in Ref. p8 



^s^2g = ^7f • (24) 

For 0.5 < Ax < 3, the decay rate is IQ-^^GeV < T^^^g < IQ-^^GeV with a 
comparable decay rate into two gluinos. This gives the saxino an approximate 
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lifetime of 10~^s < Tg^^q < 10~^s. The heavy axino and saxino have very 
short hfetimes and decay before nucleosynthesis starts. Their decay therefore 
does not affect the standard nucleosynthesis calculations. 

Finally, the relic gluino density, assuming the gluino is the LSP, should 
be small enough that it does not constitute a significant fraction of the dark 



matter halo density (see for example |2^). Refs. 0, |TT| show that taking 



into account the non-purturbative effects of gluino-gluino annihilation into 
quark- ant iquark and gluon-gluon the relic gluino density is extremely small; 
flgh"^ ~ 10^^ — 10^^^. There are however stringent limits on a gluino LSP 



coming from searches for anomalous heavy isotopes pO[ and from energetic 



neutrinos due to gluino annihilations in the sun ||31|. Both of these latter 
constraints apparently rule out an absolutely stable gluino LSP. They do 
not however constrain the case of a gluino NLSP and gravitino LSP also 
considered in this paper. 

6 Conclusions 

In this paper we have presented a solution to the fi and strong CP problems 
in the presence of a heavy gluino LSP. The model has a natural Peccei-Quinn 
symmetry which prevents the /j, term at tree level. However when the PQ 
symmetry is broken at the messenger scale M ~ 10^^ GeV the /i term is 
generated. 

The particle phenomenology of the model is quite novel. Either the gluino 
or the gravitino is the LSP. If the gravitino is the LSP, then the gluino is the 
NLSP with a lifetime on the order of one month or longer. In either case this 
heavy gluino, with mass in the range 25 - 35 GeV, can be treated as a stable 
particle with respect to experiments at high energy accelerators. 

We have studied some of the phenomenological constraints on the model. 
The most significant comes from LEP searches for the neutral Higgs boson. 
Our model is most like the no stop-quark mixing benchmark which is severely 
constrained by the data. In fact the model only survives in a narrow region 
of parameter space resulting in a light neutral Higgs with mass ~ 86 — 91 
GeV and tan/5 ~ 9 — 14. In addition the lightest stop and neutralino have 
mass ~ 100 — 122 GeV and ~ 50 — 72 GeV, respectively. Thus the model 
will soon be tested. Finally, the invisible axion resulting from PQ symmetry 
breaking is a cold dark matter candidate. 
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8 Appendix 

A The complete model 

In this section we present the complete model with all the symmetries and 
charges. The model is defined at the GUT scale by the SO(IO) invariant 
superpotential W D Wi + W2 + W3 + W4 + W^ and a non-renormalizable 
term in the Kahler potential K where 

Wi = 16310^163, 

W2 = XalOHAlOA + AxXlOi, 

W3 = Xif]iAr]i + X2f]2Ar]2 + \Xf]ir]2, (25) 

Wi = Ay3r7iri63, 

K D ^10^ + h.c. (26) 

(163,r7i,r72,'?/', V'O are all 16's, {f]i,f]2,4',ip') are 16's, (10//, 10a) are lO's, 
(X, Y) are singlets and (A) is an adjoint under SO (10). Wi, W2 and W3 were 
discussed earlier in section ^ 

The theory is invariant under U(l) PQ and R symmetries. The charges 
of the fields under these symmetries are given in Table 0. 

The R symmetry is broken by the vevs of several different fields at the 
GUT scale. However, the PQ symmetry is not broken at the GUT scale and 
prevents a /x term in the superpotential. 
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fields 


IOh 


10a 


A 


X 


Vi 


Vi 


V2 


^2 


I63 


R 


+ 1 


+ 1 


+2 


+2 


+1 


+ 1 


+ 1 


+ 1 


+3/2 


PQ 


+ 1 


-1 





+2 


-1/2 


+ 1/2 


-5/2 


+5/2 


-1/2 



fields 


^ 


^ 


v^' 


^' 


Y 


R 








-1/2 


+3 


+3/2 


PQ 








+1 


-1 






Table 1: The R and PQ charges of different fields in the complete model. 

W4 contains the field Y which gets a vev 

{Y) = Mg. (27) 

W4 also results in a bottom-r Yukawa coupling non-unification at the GUT 



scale as discussed in B.l. 



We use W5 in a standard way to split the top and bottom Yukawa cou- 
plings by giving vevs to ip and ip of order Mg- This mechanism is discussed 
in^. 

B Yukawa coupling non-unification in an SO (10) 
SUSY GUT 

In minimal SO (10), all standard model fermions in a given generation are 
contained in a single spinor (16) representation of SO(IO). The coupling of 
the form Wi results in a unified Yukawa couplings for the top and bottom 
quarks and the r lepton, i.e. Aj = A;, = A^ at Mg- It is interesting, and 
necessary for our model, to see if it is possible to relax this condition in a 
simple way. The low tan/3 fit of our SO(IO) SUSY GUT to the infrared-scale 
physical observables requires Xt — Xb splitting at the GUT scale. We are 
also interested in the possibility of splitting A^ — A^. The reason is that our 
best fits to the data come from bottom Yukawa couplings which are ~ 30% 
smaller than the r Yukawa coupling at the GUT scale. 
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B.l Xb — Xr non- unification 

Wi, W3 and W4 contain interaction terms 

miXiAr], + Ay3yi63 + AXr/2), (28) 

resulting in a heavy (16") and a massless (at the GUT scale) (16') multiplet 
given by 

16"cx(Ai(A)r7i + Ay3(F)163) (29) 

and 

16'(x{\Y3{Y)r]i-\i{A)16,). (30) 

Note, we can safely ignore the last term in Eq. ^ since the vev of X is much 
smaller than Mq- 

The massless 16' is identified with the matter multiplet containing the 
third generation quarks and leptons. Since A gets a vev in the B — L direction 
and B — L quantum numbers are different for quarks and leptons, X^ — 
At- unification will be lost, but we still have A^ — Af, unification. In fact 
this mechanism is an SO(IO) version of the one introduced in Ref. [B^ . 
Considering Eqs. H, ^ we define 

rj.B-L 

(sin eB-L)a = I " (31) 



where T^~^ is the B — L quantum number of the state a = Q,U, D, L, E, v 
in the matter multiplet and p = Xi/Xy^- 

From Wi we see that the fermions in the matter multiplet 16' obtain 
mass at the electroweak scale due to the term 



We therefore have 



(-5,16') 10^ i-SaW) . (32) 

Xb SqSd 1 1 + p^ 



A, slSe 91 + pV9' 



(33) 



^^In calculating the one loop gaugino (two loop scalar) mass contributions from W3 in 
Eq. ^ (|l^), we have ignored an order 1/p correction coming from W4. Since 6 is a small 
free parameter, this does not affect the wino, bino and scalar masses while the gluino mass 
is varied by the free parameter b. 
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where A;, and A^ are the effective bottom and r Yukawa couphngs. We thus 
find 

i<^<l. (34) 



depending on the choice of p = Ai/Ayg. For p ~ 4.2 we get \ ^ ~ 30% 



Ar 



0. 



B.2 Ai — Xh non- unification 

The splitting between At and A^, A,- Yukawa couphngs is best achieved in 
the Higgs sector using W^. Since ip and ij) get vevs of order M^, the Higgs 
doublet mass term can be written as 



( du dA d^,)\ Xx{X) Ma • (35) 

V A^(V') Xy,{Y) j \d^, J 

From the above we see that the two light Higgs doublets are dn and cos 'jdn — 
siwjd^' which are identified with the light Higgs doublets of the MSSM, 
{Hu, Hd), respectively, cos 7 is given by 

cos7 = ^^ (36) 

where 

Note, the rest of the Higgs doublets remain very heavy. From Wi we see that 

— = cos 7 (38) 

Ai 

where A^ and A;, are the effective Yukawa couplings of the top and bottom 
quarks. A hierarchy of 50 is easily achieved by choosing p' ~ 1/50. 



C T neutrino mass 

In this section we show that it is possible to get a reasonable r neutrino mass, 
in agreement with atmospheric neutrino oscillations. Let us add two SO(IO) 
singlets A^ and P to the model with R and PQ charges given in Table ^ 



fields 


N 


P 


R 


+5/2 


-1 


PQ 


+ 1/2 


-1 



Table 2: The R and PQ charges of the new fields in the neutrino sector. 



The only possible couplings for these singlets are 



(39) 



We assume that P gets a vev of order the messenger scale (10^^ GeV). The 
neutrino mass matrix is then given by 



( i^L i^R N 



( 



uit 



nit 













\ 



V -SEXnA^) ^n,{P) J 

The mass of the lightest neutrino, identified as u^, is given by 



^R 

N 



(40) 



(41) 
Thus one gets a reasonable value for rUy^ with A„j ~ 3 x 10 ^ and A^^ ~ 1- 
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Figure 1: —fi at the messenger scale, — /i, —B/fi and vB at the EWSB scale 
are plotted in this figure. 
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Figure 2: Variation of tan/3 versus d for fixed A = 10^ GeV. 
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Figure 3: MSSM Higgs masses versus d for fixed A = 10^ GeV. 
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Figure 4: A magnified variation of tfie mass of the ligfitest neutral Higgs, h 
for small d and A. The solid line shows the variation versus d for a fixed 
A = 10^ GeV while the dashed line shows the variation versus A with a fixed 
d = l. 
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Figure 5: Variation of tan/? versus A for fixed d = 1. 
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Figure 6: The Chargino and neutralino masses versus d for fixed A = 10^ 
GeV. 
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Figure 7: Squark and slepton masses versus d for fixed A = 10^ GeV. 
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Figure 8: Squark, slepton, chargino, neutralino and Higgs masses versus A 
for fixed d = \. Many of the masses are almost degenerate, therefore one 
representative is shown from each set. 
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Figure 9: The ratio of the SUSY contribution to the amphtude of 6 ^ 57 
to the SM contribution. The dashed hne shows the variation versus d for a 
fixed A = 10^ GeV while the sohd hne shows the variation versus A with a 
fixed d = 1. 
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